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The problem of two interacting particles in a quasiperiodic potential is addressed. Using analytical 
and numerical methods, we explore the spectral properties and eigenstates structure from the weak 
to the strong interaction case. More precisely, a semiclassical approach based on non commutative 
geometry techniques permits to understand the intricate structure of such a spectrum. An interac- 
tion induced localization effect is furthermore emphasized. We discuss the application of our results 
on a two-dimensional model of two particles in a uniform magnetic field with on-site interaction. 

PACS numbers: 05.45.+b, 72.15.Qm, 72.10.Bg 



I. INTRODUCTION 

The study of crystal electrons submitted to a mag- 
netic field has been-extensive since the early works of 
LandauEl and Peierlso. These studies have led to deep in- 
sights in the physics of electrons in. solids (interpretation 
of the de Haas van Alphen effectB, investigation of the 
Fermi surface...). The number of contributions on the 
subject between 1950 and 1970 reveals the importance of 
magnetic field effects. Twenty years ago, Hofstadter nu=. 
merically computed the spectrum of the Harper modelo 
and discovered its fractal structure as aJhnction of the 
normalized magnetic flux per lattice celltJ (Fig. 1). 



FIG. 1. Hofstadter's butterfly for rational values of 
a = p/q up to g = 29. For each value of the magnetic 
flux a = p/q, we generally have q bands. Near energies 
equal to ±4 and zero flux, we observe the emergence of 
Landau levels. 



The problem of a two dimensional electron on a pe- 
riodic lattice has been of special interest in solid stat, 
physics during the last fifteen years : superconductin, 



have 




and normal-metal networks^. Harper-like mod; 
been used to describe the quantum Hall effee 
ganic conductors, in Anyon supecconductivitj 
flux phases for the Hubbard modelEJ. 

If the lattice is given by the positions of the ions of 
a metal, the lattice spacing a is of the order of 1 A, so 
that even with the highest magnetic fields that can be 
produced now, namely B « 20T, we get a — j/2tt w 
0.5 x 10 -4 which is fairly small and shows that in this 
situation a "semiclassical" approximation will always be 
relevant. As a matter of fact, an effective Planck's con- 
stant denoted by 7 proportional to the applied magnetic 
field naturally appears as an adjustable variable of the 
problem. Therefore the weak magnetic field limit 7 1— > 
corresponds to the semiclassical limit h t— > 0. The cor- 
responding classical phase space at B = is the quasi- 
momentum space, namely the Brillouin zone of the cor- 
responding lattice. Topologically it is a 2-torus and the 
appearance of the magnetic field transforms it into a non 
commutative 2-torusEi 

Whenever 7 = 2np/q, (p,q 6 N) the lattice Hamil- 
tonian H recovers some periodicity and Bloch's theory 
applies. We shall see then that H can be represented 
as a self-adjoint q x q matrix whose entries are periodic 
functions of the quasimomentum components. Thus, if 
7 is close to any rational multiple of 27r, it is possible to 
compute the spectrum using semiclassical methods. 

Based on these remarks, many theoretical and math- 
ematical works were published during thc-Jast fifteen 
years using renormalization group-analysistj and pseu- 
dodifferential operators techniques!!^. On .the basis of the 
techniques of non commutative geometrytll, another ap- 
proach was developed in-order to reformulate and extend 
the semiclassical resultsEJ. The algebraic semiclassical 
approach is justified by the simplicity of its application 
and its efficiency, for example in the computation of Lan- 



and in a model- 
The comparison 



dau levels both in Harper-like mod' 
Hamiltonian on a triangular latticed 
between semiclassical formulae and exact calculations ex- 
tracted from the various spectra for 7 G 27rQ gives sur- 
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prisingly accurate agreement even for relatively large 7's 
(namely j/2n < 0.2). 

While in the above formulation of the problem of Bloch 
electrons in a magnetic field the particles are considered 
on a two-dimensional lattice, it is possible to map it ex- 
actly onto a one-dimensional lattice with quasiperiodic 
potential. The interesting property of such a lattice is 
the duality between momentum apd spatial coordinates 
pointed out by Aubry and Andrai3. This Aubry dual- 
ity results in a delocalized structure of the eigenstates 
characterized by an algebraic decay and a multifractal 
eigenspectrum. This leads toiSumasidiffusive wave packet 
spreading on such a latticeHlia. 

Recently, numbers of authors have followed a new path 
in the study of the combined effect of interaction and 
disorder. The a priori simple ppablem of two interacting 
particles in a random potentially has indeed revealed an 
unsuspcctly large interaction induced derealization ef- 
fect. However, the opposite effect has been discovered in 
the case of two particles in a quasiperiodic potential. In 
this case, the interaction leads to the emergence of a pure- 
point component out of the spectrum of the non inter- 
acting problem. These facts have been firmly establ ished 
by overconvincing numerical and analytical resultsr 1 lr 2 l. 
It is one of the purposes of this paper to again express 
these arguments in more details. 

We shall present in this work analytical and numer- 
ical results derived from the two-particle Harper prob- 
lem with on-site interaction on a one-dimensional lattice. 
More precisely we devote the second section to the pre- 
sentation of the algebraic semiclassical approach on the 
non interacting problem U — 0. The corresponding spec- 
trum is somehow an intricate superposition of two Hof- 
stadter butterflies. The aim of section 3 is to study the 
small interaction regime where usual perturbation theory 
can be applied. The evolution of the spectrum as a func- 
tion of the strength of the interaction will be presented. 
After building the analytical framework in section 4, we 
apply it to the computation of the levels in the strong 
interaction regime. We show that for very large U, the 
spectrum is divided into two parts : one corresponding 
to the non interacting case and the second one, look- 
ing like a Mathieu spectrum corresponding to localized 
states strongly influenced by the interaction. Based upon 
Aubry's dualityliZl, it can be proved that all the wave func- 
tions are localized in this regime as far as the Mathieu 
part of the spectrum is concerned. Finally, we discuss in 
section 5 the problem of two interacting particles on a 
two-dimensional lattice submitted to a magnetic flux. 

II. NON INTERACTING MODEL 

In his 1930's study of the electronic diamagnetism of 
metals, Landau computed the energy spectrum of a free 
electron subject to a uniform magnetic fieldEl. If B is 
uniform and parallel to one axis, for example axis 3, the 



kinetic energy is written as : 

with = (P,, - q e A^) /%, fi= 1,2 and A = (A 1 ,A 2 ) is 
the vector potential satisfying curl(yl) = B, q c is the elec- 
tron charge. Moreover, the quasimomenta K\, K% satisfy 

[K 1 ,K 2 ]=iq e B/h (2) 

Let us notice that this commutation rule becomes canoni- 
cal when replacing Ti by q c B/h. This new effective Planck 
constant (divided by 2ir) is proportional to the magnetic 
field B and behaves as a varying physical parameter, 
quite naturally. 

The spectrum of Hi is E n = EQh c ffU}(2is + 1) with 
Eq = h 2 /2m c , h c ff — q c B/h and oj = 1. Therefore : 

E v = hu c (v +1/2) (3) 

where oj c = q c B/m c is the cyclotronic frequency and v is 
the Landau quantum number. 

When B = 0, the electron energy E(k) for each 
conduction band is given by Bloch's theory, where the 
quasimomentum components k = (fci,fc2) are defined 
modulo the reciprocal lattice such that for a simple 
square lattice in the tight-binding approximation E(k) = 
2Eq [cos(kiai) + cos(fc2a2)] where a M is the vector of the 
Bravais lattice in the //-direction. The charge carriers en- 
ergy is calculated by expanding E{k) near its extremum, 
denoted by k c , namely : 

E(k) = E(k c ) + ti 2 (M- 1 ).. kikj/2 + O (\k\ 3 ) (4) 

where M stands for the effective mass matrix such that 
Af- 1 = D 2 E{k c )/h 2 . 

Thus Landau theory leads to a substitution ki ■ a; 1 — > 
Ki = i (P — q c A) ■ cii when an external magnetic field is 
applied. We have the following commutation rule : 

[Ki,Kj] — iq c Ba,iaj /h — 2iTT(f>ij/(/)Q = 2iira = i"f (5) 

where 0o — h/q c is the flux quantum, </>y is the magnetic 
flux through the cell generated by (at, dj) and a — <pij/4>o 
is the normalized magnetic flux. For a crystal with peri- 
odic spacing, the PeierlSr<pperator V(k) is represented by 
an effective Hamiltoniara, namely : 

7>(£0 = $> m (aK m - fc , m e Z 2 (6) 

m 

where h m (a) are smooth functions of a. Thus : 

H e v{k x ,k 2 )= h m j mk (7) 

1 a 

If several bands intersect the Fermi level, the interband 
coupling due to the magnetic field is neglected and there- 
fore : 
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H cS = 2t (cos K x + cos K 2 ) (8) 

where t is physically interpreted as a transfer term corre- 
sponding to the required energy for an electron to jump 
from one site to another (nearest neighbour) of the lat- 
tice. 

For a wave function ip(ni,n 2 ) defined on the two- 
dimensional lattice £ 2 (Z 2 ), the magnetic field effect can 
be seen through the magnetic translation operators such 
that : 

(Uiip) (ni, n 2 ) — e J(«i-i,« 2 ) tfj(n 1 — l,n 2 ) 

{U 2 ip) (ni,n 2 ) = e n J («i.^-u V(«i 5 «2-1) (9) 

in an appropriate gauge we get : 

(Kiip) (ni,n 2 ) = - !> n 2) 

(W 2 V) (m.na) = e-^(ni,n 2 - 1) (10) 

Because of the presence of the uniform magnetic field, 
the magnetic translation operators no longer commute, 
namely in that case 

U X U 2 = & ll U 2 U x (11) 

where 7 is the normalized magnetic flux per lattice-cell 
defined by 7 = 2ira — 2it<f>/(j)Q, <fi being the flux per unit 
cell and 0o = h/q e the flux quantum. 

If we set 1t\ — exp yilCxj , U 2 = cxp (iK 2 ^j using the 

commutation rule (11), we obtain 

[k u k 2 ] = = 2in ± = = i7 (12) 

n 0o 

which corresponds to (5) in the particular case i = 1 and 

3 = 2- n 

Following Harpcru, the eigenvalue equation is written 

E [tp(ni + a, n 2 ) + tp(m - a, n 2 ) + 
+\e^ Bnia/n ^(n 1 ,n 2 + a) + Xe-^ Bnia/h ^( ni ,n 2 - a) 

= 28i/j(n 1 ,n 2 ) (13) 

A represents the strength of the quasiperiodic potential. 

Let us assume plane-wave behaviour in one direction, 
i.e. we set if>(ni,n 2 ) = J dj3e l ^ n2 4>{n\) since the coeffi- 
cients in the previous equation only involve n\ : 

^(n 1 ,n 2 ) = e^cj ) (n 1 ) 

and the eigenequation becomes : 

4>(nx + 1) + (j){nx - 1) + 2A cos (2nani + f3) 0(m) 

= £4>{n 1 ) (14) 

where we included the additive energy due to the motion 
in the field direction in the eigenvalue £ and where we 
changed the origin of n\. 



It is possible to characterize the properties of eigen- 
functions from (14) by looking at a special regime, 
namely A C 1. Therefore, the hopping term is domi- 
nant and we can treat the quasiperiodic potential part 
of the eigenvalue equation as a perturbation. It is then 
easy to see that the solutions are given for A = by Bloch 
waves 4>k(n) = exp(ifcn) with an energy E — 2cosfc. For 
< A <C 1, the perturbation theory allows us to perform 
an expansion of eigenvalues and eigenstates in A such that 

E{k) = 2cosfc + ^A m e m (fc) 

m 

Mn) = e lkn (l + A m / m ( 7 n + 0)j 

= e lkn u m ( in + (3) (15) 

Evaluating the first and second order perturbation theory 
contributions and replacing the expressions (15) in (14) 
leads to : 

A (u rn+ i + u m _i) + 2 cos( 7 m + k)u m = E(k)u m (16) 

The previous equation is known as the "almost Mathieu" 
eigenvalue equation and the argument above is the Aubry 
dualityO between momentum and coordinate represen- 
tations. As far as spectral properties are concerned one 
can be easily convinced that dealing with Bloch states 
means that the states are extended. Thanks to this du- 
ality, A <-> 1/A between (14) and (16), it is quite natural 
to get localized states for the almost Mathieu Hamilto- 
nian at small A's. More precisely, it has been proved that 
the spectrum of the almost Mathieu Hamiltpnian is pure- 
point at small A's and for almost all /3'scil Conversely 
if A ^> 1, the almost Mathieu Hamillanian has purely 
continuous spectrum for almost all /3'sEj. 

Setting t = 1 in formula (8) and using the magnetic 
translation operators U\ and IA 2 defined on the two- 
dimensional square lattice by (10), the previous Harper 
equation can be written as the action of an effective 
Hamiltonian such that : 

H eS = Ux + Wf 1 + U 2 + U 2 Y (17) 

In order to study the two interacting particles model on 
a quasiperiodic lattice we transform the previous eigen- 
value equation (14) into (A = 1) 

[2 cos (7711 + Pi ) + 2 cos ( 7 n 2 + (3 2 ) + U5 ni ,„ 2 ] cj) ni , n2 

+0ni+l,n 2 + 0ni-l,n 2 + 0ni,n 2 -l-l + 0m ,713—1 

= Ecj) num (18) 

where /3i j2 are related to the quasimomentum compo- 
nents of the non interacting case. In the following we 
shall consider /3i. 2 = (3. Here we chose the form of on-site 
interaction which only influences the symmetric configu- 
rations while the antisymmetric ones remain not affected 
by U. Due to that, we shall only discuss symmetric con- 
figurations in the following. 
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In the most simple case of non interacting particles 
(U = 0), the spectrum can be computed as before and is 
shown in Fig. 2. 



FIG. 2. Spectrum of the two-particle Harper problem 
with U = obtained for rational values of a — p/q up to 
q = 19. 

As we pointed out before, 7 = 2na appears in our 
problem as an effective Planck constant since the mag- 
netic translation operators U\ and U2 obey canonical 
commutation rules in 7. Therefore, we study the semi- 
classical limit by letting 7^0. It is also possible to per- 
form calculations near a rational value of the magnetic 
flux, namely 7' = 7 — 2np/q 1 — > 0. The efficiency and 
the accuracy of our calculations allow to explain some 
features of the corresponding spectra. 
When 7 = 0, corresponding to B = 0, we recover the 
band function E(k), where k = (fcijfe). To study the 
Landau levels, we expand the classical symbol of the 
Hamiltonian around an extremum of the band function 
denoted by k c : 



n{k) = n(k c ) + ^d,H(k c )k^ + 



(19) 



The quantization of TL(k) consists in replacing the mag- 
netic translation operators bjO: 



Uj = exp (i(k cj + yfyKj)) , j = 1, 2 



(20) 



where k C j are the bottom well coordinates and Kj are 
operators satisfying Heisenberg's commutation relations 
\K\,K-i\ = i. The quantized of Tt, denoted by H, is 
written as : 



i(m-k c -\-^/ : ym-K) 



H = ^/i(m,7)e 



with to • K = ui\K\ + TO2i^2- In the weak field limit, one 
formally expands H in powers of y/j : 

H = {Km, 0)e im k ' + i^h(m, 0)e m k 'm ■ K 



dh 



(771, 0)e imk " 



_c?7 
+0( 7 3/2 ) 

which we rewrite as 



-h{m,{))e lm k "{m- Kf 



(21) 



on. 



H = H{k c , 0) + 7 ( — (kc, 0) - ~d„d u H(k c , 0)K^K v 



+0( 7 3/2 ) 



(22) 



The dH/dj-term takes into account a possible explicit 7- 
dependence of the classical Hamiltonian whereas d^d v TL 
represents the inverse effective mass matrix due to the 
band function curvature. By a unitary transformation, 
the quadratic term can be written as ui (K,\ + JC 2 ,) /2 
where cu is related to the determinant of the Hessian 
matrix d ll d 1 /H(k c ,Q). We recognize here the harmonic 
oscillator Hamiltonian. For this reason, the energy levels 
denoted by E v are called "Landau levels" and are equal, 
to that order in 7, to lo[v + 1/2) leading to : 



£,.(-. ) = HUv.U) + -(2;/ I ) ( det ] -D' 2 H(k,,0) 



( dH(k c ,0) 



0( 7 



1/2 



(23) 



The formula (23) has been checked numerically on sev- 
eral models. To illustrate it, let us consider the two- 
particle Harper Hamiltonian on the square lattice (18) 
near the maximum k c = (0, 0) of the band function. Us- 
ing (20) and (22) the quantized Hamiltonian is then ex- 
pressed as an expansion in powers of 7 : 

H = 8 - 7 ((K^r + {K^f + [K^f + 

+ ^((^ 1) ) 4 + (^ 1) ) 4 + (^ 2) ) 4 + (^ 2) ) 4 ) 



+0( 7 3 



(24) 



where the K^ 1,2 ^ are quasimomenta for particle 1 and 2 
respectively. Finally it gives the Landau levels : 



E„ 



8-27(7/1 +7/ 2 + l) 

+7 2 [(2i/i + I) 2 + (2t/ 2 + l) 2 + 2] /16 



+0( 7 3 



(25) 



where v\ and 7/2 are the Landau quantum numbers asso- 
ciated with particle 1 and 2 respectively. To check the 
accuracy of this formula, we compared it to the datas 
extracted from the numerical spectrum obtained by ex- 
act diagonalization. Fig. 3 shows the accuracy of such 
a semiclassical expansion in the description of the spec- 
trum of the two-particle Harper model when 7^0. 
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FIG. 3. Comparison between semiclassical calculations 
(25) (full curves) and exact numerical spectrum (points) 
for Landau sublevels in the two-particle Harper model 
on the square lattice when U = 0. Datas are extracted 
in the region of energies corresponding to the maximum 
(0, 0) of the band function. 



III. WEAK INTERACTION REGIME 

We present here a simple perturbative treatment that 
enables to implement the already presented results for 
the weakly interacting case. The first-order contribu- 
tion allows to understand the splitting of Landau bands 
at sufficiently weak interaction, and describes it qual- 
itatively well. It moreover enlightens the mechanism 
through which interaction affects the system. Using the 
representation defined by (20), we write the unperturbed 
Hamiltonian as : 

ffeff = 2 cos(V7^i) + 2 cos(V7^2) (26) 

In the semiclassical limit 7 h wc expand (26) in 
a power series around a minimum of potential qN — 
7r/^/7 + 2ttN/^ 1 N e Z. Keeping only terms up to 



the second order in 7 we end up with a harmonic oscilla- 
tor. In this approximation and in the continuous case the 
one-particle wave functions of the unperturbed Hamilto- 
nian are therefore given by : 

My) = M^=) cx P (-f^) /yW^ (27) 

Here, H v (x) is a Hermitc polynomial, the index v refers 
to the Landau level, y = x— qjy in term of the minimum of 
potential qN around which the harmonic approximation 
has been performed, and x is the spatial coordinate. This 
expression is of course valid, provided 7 and \x — qw\ <C 1, 
i.e. in the small magnetic field regime, and not too far 
away from a potential minimum. Extending our expan- 
sion to higher powers in 7 would allow us to increase the 
range of validity of this expression. We could indeed write 
the exact normalized wave functions in an expansion in 
7 as 

MV) = cxp (-f^) (c H„(j=) +7Ci^ 1) (J=) + ...) 

(28) 

For the purpose of discretization, wc introduce a con- 
tinuous variable £ € R labeling the well, and a discrete 
one I G Z numbering the sites. Then y = £ — Z^/7 
since in the chosen representation, the intersite spacing 
is a = y/j. The set {<fiv} builds a quasiorthogonal basis 
in the sense that for £ 7^ due to the Gaussian envelop 
of the states we have : 

5>"(* " 1>/T)MZ' - l V^) = 0(exp(-l/7))<W (29) 

These functions are periodic in £ with period V-^/7- In 
the semiclassical limit the norm of </?„ is : 

00 

\\M\ 2 = E M£-V7)l 2 

I — — 00 

= 1/V7 J dy\Mv)\ 2 

= VV7 (30) 

Consequently, to get normalized one-particle wave func- 
tions on the discrete lattice £(Z) we must multiply the 
tp's by a factor 7 1 / 4 . We thus can write the symmetrized 
two-particle unperturbed wave functions as : 

±M€ - iWM? - l 'i)) (1 - <W(! - VV2)) (31) 

We are now able to compute the first-order correction 
for the energy. Because of the exponentially localized 
character of (28), two particles located on different wells 
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have only an exponentially small overlap, and as a con- 
sequence do practically not interact. Therefore the first- 
order interaction induced correction to the energy is non 
zero only for symmetric wave functions with £ = We 
have : 

aeW = tf£ $>«'(z, z ')) 2 <5(£ - £' + (*' - OVt) 

i 2' 

= f/ %£' y dy (<pn(y)<p v (y)) 2 (2- <J M ,„) + 
+0(cxp(-l/ 7 )) (32) 



for states with Landau quantum numbers (0,0) and (0,1). 
This result shows that the interaction primarily acts on 
two-particle states with high double-site occupancy. In 
what follows we shall call such states "pair states" . States 
for which the particles are located around different poten- 
tial minima practically do not feel the interaction. There- 
fore, switching on the interaction does not modify most 
of the spectrum as can be seen on Fig. 4. 

From (25) and (34) and for small enough interaction, 
the shifted part of the spectrum is given by : 

E„ u u 2 (7) ~ 8 + U x [^- - 2 7 (n +2^2 + 1) 

+ 7 2 [{2vx + if + (2u 2 + l) 2 + 2] /16 (35) 

The amazing agreement between the numerically com- 
puted spectrum obtained by exact Lanczos diagonaliza- 
tion and (35) is shown in Fig. 5 where U = 0.4. It is 
a confirmation of our reasoning : pair states form the 
shifted part of the spectrum. Because these states are 
much fewer than states where particles are located in 
different wells, the shifted spectrum is much less dense. 
In this sense the interaction splits the butterfly into two 
parts. One of them is practically not affected by the in- 
teraction and corresponds to the states where particles 
are far from each other. The second one is shifted and 
relays to the situation where particles form pair states. 
Here, the interaction results in a global shift of the spec- 
trum. In this way, new states appear in the initial gaps of 
the-non interacting spectrum (see Figs. 4,6 and Fig. 1(b) 
ir£d). Direct analysis of eigenstates shows tiiat the cor- 
responding states are exponentially localizedEil. We shall 
come back to this point later on for the case of strong 
interaction. 



FIG. 4. Spectrum of the two-particle Harper model 
with on-site interaction at U — 0.4 up to q = 23. 

From (28), the dominant term in the last integral is of 
order O(-Jj) so that we finally have 

A£ (1) - US u >y/j (33) 

The numerical factor can be estimated from the harmonic 
approximation (27) which leads to : 



A^ 15 = U6 U > J^ = U5 U , (34) 



G 



butterfly (see Fig. 4 where U = 0.4), then evolving to a 
shifted Mathieu spectrum as the interaction grows bigger 
and bigger (see Figs. 6, 7 and 8 where U = 5, 10 and 20 
respectively). 



FIG. 5. Comparison between semiclassical calculations 
extended by perturbation theory (35) (full curves) and 
exact numerical spectrum (points) for the two-particle 
Harper model with on-site interaction at U = 0.4. 



IV. STRONG INTERACTION REGIME 

The strongly interacting regime needs a special treat- 
ment quite analogous to the one presented in section 2. 
As we will see, Schur's complement formula can be suc- 
cessfully applied to construct an effective Hamiltonian. 
The latter is then expanded in a power series in 7 to de- 
liver highly accurate formulae. From the weakly interact- 
ing regime we learned that particles located on different 
potential minima do not feci each other : for such pairs, 
the interaction is suppressed by an exponentially small 
term of order 0(t/exp(— 1/7)). Therefore, this picture 
remains valid even for large t/'s, the relevant parameter 
being the magnetic flux. Pair states on the other hand 
undergo an energy increase of order AE w U. Therefore 
when the strength of the interaction U > increases, 
one part of the spectrum is almost not affected. Another 
spectral structure appears, initially looking like a shifted 



FIG. 6. Spectrum of the two-particle Harper model in 
the intermediate regime U = 5 up to q — 23. 
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FIG. 7. Spectrum of the two-particle Harper model in 
the strongly interacting regime U — 10 up to q = 23. 



FIG. 8. Spectrum of the two-particle Harper model in 
the strongly interacting regime U = 20 up to q = 23. 

In this section, we present an analytical approach that 
allows to understand completely the mechanism driving 
this evolution of the spectrum. Further details like the 
splitting of the Landau band v\ = 0, = 1 will also be 
computed, even though the physics is there less transpar- 
ent (see Fig. 7). We shall concentrate our semiclassical 
calculation near the band function maximum k c = (0, 0) 
corresponding to the energy z w U + 4 in the spectrum. 
The two-particle Hamiltonian can be expressed in the 
following way : 

#tip - E m ,„ I 2 cos h m + P) + 2 cos {in + P)\ 
\m <g> n)(m <g> n| + UJ2 m \ m ^rn}(m<^) m\ 

(36) 

+ Em/n |m (g) n) [(m (g) n + 1| + (m (g) n - 1| 
+ (m + 1 O n\ + (to - 1 <g> n\] 

The strategy is based on the so-called Schur complement 
formula. Our Hamiltonian ifrip is a self-adjoint operator 
acting on a Hilbcrt space that can be decomposed as 
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H = V ® Q. Let P and Q be the orthogonal projections 
on each subspace of that decomposition, namely : 

P = J2 m \ m ®fn){fn® m\ 

Q = I - P = J2 m jkn \m (g) n) (m (g) n\ 

In other words, P is the eigenprojection on pair states 
and Q is its orthogonal. If z is an eigenvalue of ifrip and 
does not belong to the spectrum of QHtipQ then it is 
also an eigenvalue of the following effective Hamiltonian 

Hf IP (z) = PH T1P P + PH T1P Q J tt -QHt 



z — QHtipQ 



2 TIP 



P 



(37) 



When U is large the dominant term in the effective 
Hamiltonian given by the Schur complement formula (37) 
corresponds to the pair states. The semiclassical ap- 
proach we introduced in section 2 remains valid so that 
Hf IP (z) = Hf lP (z a + 701 + 7 2 z 2 + 0( 7 3 )). The implicit 
equation to be solved is then : 

H T %(z) = zo+7^+7 2 22 + 0( 7 3 ) (38) 

with 

<W = 4°ipW + iH&iz) + 7 2 4ip(^) + 0(7 3 ) 

(39) 

The expansion of the dominant term reads : 
PHnpP = U + 4 cos(^/jK 2 ) = U + 4 - 2-yK% 
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(40) 



K\ + 0( 7 d ) 



and if we consider U large, z is large too so that : 



1 



z — QHtipQ 



1 QHtipQ QHtipQQHtipQ 



+0 (z- 4 ) 



(41) 



Expressing the different contributions in Schur's formula 
and expanding in powers of 7 lead to : 



H (0) 
-"tip 



„ TT 8 32 176 , _ 4N 
4 + U+ - z + ^ + — + °( z ) 



H T %(z) = 



-2(z 3 + 8z + 64) 



K. 



z 3 + 8z + 256 



A' 



z 2 + 4z + 34 g2 
2 ^z 3 + 8z + 64 1 

z 2 + Az + 70 4 ' 
z 3 + 8z + 256 1 

,z + 8 



+ 16 



(z + ; 



z 3 (z 3 + 8z + 64) 



(AT? + Iff 



(42) 



Finally, we have to solve (38) to get the coefficients zq, 
z\ and Z2- The corresponding equations for those coeffi- 
cients are at most of degree four. We shall give here the 
equation that z has to satisfy at the order O [z~ ) 



4 + U + 



32 176 



(43) 



In a very similar way used for the computation of z , 
the analytical expressions of z\ and z-i can be derived 
from (38), (42) and (43). The good agreement with the 
exact numerical spectrum can be seen on Fig. 9 for 
U = 50. Here the numerical values for the sublevels 
are : for i/ 1)2 = 0, 2(7) = 54.1597 - 0.2826 7 + 0.03567 2 , 
for ^1,2 = (0,1), 2(7) = 54.1597 - 0.84807 + 0.20847 2 , 
for z/ M = (1,1), z{i) = 54.1597- 1.41337 + 0.5539 7 2 . 
A similar computation can be done near the band func- 
tion minimum k c = (n, n) corresponding to the energy 
zksU-4. 



FIG. 9. Comparison between semiclassical calcula- 
tions (full curves, see text) and exact numerical spectrum 
(points) for levels in the two-particle Harper model for 
U = 50. 

The structure of the pair states for U » 1 can be 
understood in the following way : the diagonal term cor- 
responding to the energy of particles located on the same 
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site is 4A cos(7n + 0) + U. The transition amplitude on 
the diagonal ni,2 = n is given by the amplitude of the 
hopping via virtual states with n\ — n% = ±1 and energy 
denominator 1/U. There are two such paths so that the 
effective amplitude is V c s — 2/U. The same expression 
can be derived by the Schur formalism (see Sec. IV). 
After dividing the Hamiltonian by V e fi we arrive to the 
eigenfunctions equation in the form of Harper (14) with 
A replaced by A G fr — U ^> 1. Since A G fr ^> 1 when U ^> 1, 
the pair states are always within the localized phase of 
the Harper equation showing exponential localization. In 
Fig. 10, we show a typical eigenstate of the Mathieu part 
of the spectrum for U — 50 and 7/2-77 = 34/55. The fact 
that it is localized confirms the pure-point character of 
the corresponding spectrum. 

Above we showed that in the case of strong interac- 
tion, we have A e g 3> 1. This explains the appearance of 
a pure-point component in the spectrum. However, we 
think that this pure-point component will even appear 
for small values of the interaction. Our argument is the 
following : without interaction, the system obeys Aubry's 
duality while the presence of the interaction introduces 
Aubry's duality breaking. Indeed, from (18) it is easy to 
see that the interaction acts in the coordinate space and 
the symmetry with momentum space disappears when 
U 7^ 0. Formally, this argument is not sufficient to prove 
the existence of pure-point spectrum at arbitrary small 
U. Howcvc£.lhe ensemble of numerical datas we have 
here and inEj'tll confirms this conjecture. 




the case U = (25). To derive this splitting using semi- 
classical analysis, we again apply the Schur complement 
formula. Dealing with the unshifted butterfly leads us to 
consider as the dominant term QH^\pQ such that (37) 
becomes : 

Hf w {z) = QH T1P Q + QH T1P P -L -PH T1P Q 

(44) 

Applying the same scheme as before produces an ad- 
ditional shift from the unperturbed energy given in first 
order in 7 by : 

1^(7)1=4^ (45) 

This shift is valid for the second Landau sublevel (17 = 
0, vi = 1), its accuracy is shown in Fig. 11 and the two 
splitted subbands are given by : E(j) = 8 — 4.I6667 and 
-£(7) = 8 — 47 up to order 1 in 7. 



FIG. 10. Semilogplot of W = |0(U')| 2 for a localized 
state (E=50.25 -30 <\nW < -1). 

When U is large, the unshifted part of the spectrum 
looks very much like the spectrum at U = 0. The main 
difference can be found by looking carefully at the Lan- 
dau levels (see Fig. 11). The reminiscence of the exis- 
tence of the interaction is seen through the appearance 
of a splitting of Landau sublevels. This splitting only ex- 
ists when Landau quantum numbers are different v\ 7^ v<i 
and the two particles are located in the same well. Such 
a behaviour is illustrated by Fig. 11. The other sublevels 
are described by the semiclassical formulae obtained in 



FIG. 11. Semiclassical calculations (45) (full curve) 
and exact numerical spectrum (points) for the splitting 
of the V\ = 0, V2 = 1 Landau sublevel in the two-particle 
Harper model for U = 20. 
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V. TWO INTERACTING PARTICLES ON A 
TWO-DIMENSIONAL LATTICE 

Even though the studied model was derived from a 
model of two-dimensional electrons, its effective dimen- 
sion is 1 : as we already pointed out, (18) was derived 
assuming that the particle propagate as plane-wave in 
one direction. This assumption, though reasonable in 
the one-particle model, could be violated by interaction 
induced quantum interferences in the two-particle case. 
Therefore the question of the survival of interaction in- 
duced localization effect for two interacting particles in 
two dimensions remains an open problem. In this section 
we would like to discuss briefly this situation. For two 
interacting particles moving on a two-dimensional square 
lattice submitted to a uniform magnetic flux, the eigen- 
value equation reads : 

c YXi + 1, J/1, £2, J/2 ~ C TX 1 —l,y 1 ,X 2 ,y2 

Yx 1 ,y 1 ,X2 + l,y2 ^ c rii,!/i,i2-l,»2 
+V'a;i,yi,2;2,y2 + 1 ~h ' l Px 1 ,y 1 ,x 2 ,y2 — 1 
~^U5 XlyX2 5 yii y 2 ll>x 1 ,y 1 ,x 2 ,y2 — E'4>x 1 ,y 1 ,X2,y2 (46) 

where (21,2, 2/1,2) are integers denoting the positions on 
the square lattice and U is the on-site interparticle in- 
teraction. For U = 0, the previous equation can be re- 
duced to the one-dimensional Harper equation we dis- 
cussed above (18). With interaction, the same equation 
(18) can be obtained in the ansatz of plane waves prop- 
agating in one direction with renormalized interaction 
Uc3. While this plane wave approximation is a standard 
approach for the one-particle Harper problem, it has to 
be handled with care in the interacting case. Indeed this 
plane wave ansatz breaks the symmetry of the original 
problem (46). This symmetry can be seen in the limit of 
strong interaction !7> 1. In this case, there should be 
two energy bands : one corresponding to the pair states 
when particles are located on the same site with energy 
E w U and the other with E w 1 for the states in which 
the two particles avoid each other. In the higher energy 
band, the eigenvalue equation for the pair states up to 
the terms of order 1/U has the form : 

2 

— (e 2ny <j) x+hy + e~ 2l ~< V 4> x -l, y + <j>x,y + l + <f>x,y-l) 

+U<t> x . y = E<p x , y (47) 

Here the term 2/U represents the transition amplitude 
for pair states. Its derivation is similar to the case of 
two interacting particles in the one-dimensional Harper 
model. Indeed if one keeps x\ — X2 then the hopping 
term is given by V e g = 2/U because there are two paths 
with virtual energy U (2/1,2 - * 2/1,2 + 1) which contribute 
to the hopping term in the y- direction. Similarly the 
hopping in the x-direction is V c s — 2e ±2l ~ l /U. 

This representation shows that the symmetry between 
the two directions or the Aubry duality is not broken by 



the interaction. The main reason is that the symmetry of 
the interaction is invariant under rotations on the square 
lattice. In the limit of large U, this property can be seen 
through equation (47). However the symmetry ( Aubry 's 
duality) should also be preserved for small interaction. 
Due to that, we expect that similarly to the Harper model 
with A = 1, the interaction will not generate pure-point 
component in the spectrum. However this conjecture has 
to be directly checked in further analytical and numerical 
studies. 



VI. CONCLUSIONS 

In this paper we have emphasized a localizing effect 
due to the combined action of an on-site interaction and 
a quasineridodic potential. Unlike in the random poten- 
tial caseo, extended unperturbed states are localized by 
the interaction, and this localization occurs at arbitrar- 
ily small attractive/repulsive interaction. We success- 
fully identified the mechanism responsible for this effect 
as a decoupling of a Mathieu-like model from the original 
two-particle Harper model in the limit of large interac- 
tion. Our conjecture is that a similar mechanism will also 
work for small U due to an interaction induced break- 
ing of Aubry's duality. This breaking happens in one- 
dimensional incommensurate models, however in two- 
dimensional magnetic models, we expect that the inter- 
action will not break the duality and that a pure-point 
component in the spectrum will not arise. Further veri- 
fications of these conjectures are required. 
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